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ABSTRACT 

The origin of the luminosity dependence of the equivalent width (EW) of broad emis- 
sion lines in AGN (the Baldwin effect) is not firmly established yet. We explore this 
question for the broad C iv A 1549 line using the Boroson & Green sample of the 87 
z < 0.5 Bright Quasar Survey (BQS) quasars. Useful UV spectra of the Civ region 
are available for 81 of the objects, which are used to explore the dependence of the 
C iv EW on various emission properties. We confirm earlier results on the strong corre- 
lations of the C iv EW with some of the emission parameters which define the Boroson 
& Green Eigenvector 1, and with the optical to X-ray slope a ox . In addition, we find 
a strong correlation of the Civ EW with the relative accretion rate, L/LEdd- Since 
L/L-Edd drives some of the Eigenvector 1 correlations, it may be the primary physical 
parameter which drives the Baldwin effect for C iv. 
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1 INTRODUCTION 

The inverse correlation of the broad emission line equivalent 
width (EW) with luminosity in AGN, discovered by Bald- 
win (1977) for the Civ A1549 line, was intensively explored 
over the past 20 years (see a comprehensive review by Os- 
mer & Shields 1999; and more recent studies by Wilkes et 
al. 1999; Green et al. 2001; Croom et al. 2002; Dietrich et 
al. 2002; Kuraszkiewicz et al. 2002; and Shang et al. 2003). 
The physical origin for this effect is not clearly established 
yet, but a plausible explanation is softening of the ioniz- 
ing continuum shape, and increasing gas metalicity with in- 
creasing luminosity (Korista, Baldwin & Ferland 1998). The 
softening of the ionizing continuum shape with luminosity is 
supported by: 1. some observed correlations of the emission 
line strength with the ionizing spectral shape (Green 1998; 
Wang, Lu & Zhou 1998), 2. predictions of simple accretion 
disk models for the ionizing spectral shape dependence on 
luminosity (Netzer, Laor & Gondhalekar 1992), and 3. the 
dependence of the slope of the Baldwin effect on the ioniza- 
tion potential of the emitting ions (Espey & Andreadis 1999; 
Green et al. 2001; Dietrich et al. 2002; Kuraszkiewicz et al. 
2002). Non isotropic continuum emission (Netzer 1987), and 
the intrinsic Baldwin effect (Pogge & Peterson 1992; also re- 
view in Osmer & Shields 1999, §5) can produce some of the 
observed scatter in the Baldwin effect for the various lines. 

Significant correlations also exist among various opti- 
cal emission line properties in AGN, including the [O III] 
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and Fe II strength, and the H/3 FWHM, which form part of 
Eigenvector 1 (EV1) in the comprehensive study of Boro- 
son & Green (1992, hereafter BG92). Boroson, Persson & 
Oke (1985) speculated that the underlying physical param- 
eter which drives these correlation is the relative accretion 
rate, L/L^dd, a speculation which received strong support 
through recent observations (see review in Laor 2000). Fur- 
thermore, Wills et al. (1999) found that various UV emission 
properties, including the C iv EW, are correlated with the 
optical EV1 parameters, which suggests that the C iv EW 
may also be largely driven by L/L^dd- The Baldwin effect 
may then just be a secondary correlation induced by the ten- 
dency of more luminous AGN to have a higher L/L^dd- The 
possible role of L/L-^dd in the Baldwin effect was already 
mentioned by Brotherton & Francis (1999). Later, Wilkes 
et al. (1999) noted that narrow line AGN are outliers to the 
Baldwin effect, but the meaning of this result was not inter- 
preted at that time (see §4 here). Most recently Shang et al. 
(2003) suggested an indirect indication, based on spectral 
principle component analysis, that L/L-^dd may contribute 
to the scatter in the Baldwin effect. 

Recent studies established that reasonably accurate es- 
timates of the black hole mass (A/bh) can be obtained 
in AGN based on the continuum luminosity and H/3 
FWHM (M BH oc L 1/2 (H/3 FWHM) 2 , and thus L/L E d d oc 
L 1/2 (H/3 FWHM) -2 , e.g. Laor 1998). This opens up the pos- 
sibility to test directly whether the Baldwin effect is driven 
by L/LEdd, which is the main point of this paper. The data 
set and measurement procedure are described in §2, the cor- 
relation analysis is presented in §3, and the main results are 
summarized in 54. 
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2 THE MEASUREMENTS 

For the purpose of this analysis we use the BG92 sam- 
ple which includes the 87 z < 0.5 AGN from the Bright 
Quasars Survey (BQS; Schmidt & Green 1983). This sam- 
ple extends in luminosity from Seyfert galaxies with vL v — 
3.3 x 10 43 erg s _1 (calculated using the continuum fluxes in 
Neugebauer et al. 1987, assuming Ho = 80 km s _1 Mpc -1 , 
fio = 1-0), to luminous quasars at vL v = 1.4 x 10 46 erg s~\ 
where vL v is calculated at 3000 A. This is a complete 
and well defined sample, selected based on (blue) color and 
(point like) morphology, independently of the emission line 
strengths. It is also the most thoroughly explored sample 
of AGN, with a wealth of high quality data at most wave 
bands. The EV1 correlations were established by BG92 using 
this sample, and the ability to obtain reasonably accurate 
estimates of Mbh (and thus -L/I/Bdd) was demonstrated for 
this sample (Laor 1998), thus making it an optimal sample 
to explore the possible origin of the Baldwin effect. 

Archival UV spectra of the C IV region are available for 
85 of the 87 BG92 objects. The HST archives contain UV 
spectra of 47 of the objects, which were obtained by the 
Faint Object Spectrograph (FOS); the UV spectra for the 
remaining 38 objects with no HST spectra were obtained 
from the IUE archives (see Table 1). An average spectrum, 
weighted by the S/N ratio, was calculated when more than 
one archival spectrum was available. This averaging should 
decrease the scatter in the Baldwin effect induced by vari- 
ability, as demonstrated in the intrinsic Baldwin effect. 

Three of the archival spectra (PG 0934+013, 
PG 1004+130 and PG 1448+273) did not have a suffi- 
cient S/N to measure the Civ EW, and in one object 
(PG 1700+518) C IV is heavily absorbed (e.g. Laor & Brandt 
2002), leaving a sample of 81 objects with a measurable 
C IV EW. The C IV region of PG 0043+039 is also rather 
heavily absorbed, and its measured C IV EW is probably just 
a lower limit. The continuum flux uncertainty of the HST 
spectra, which were obtained without the flux density uncer- 
tainty from the HST archives, was estimated using the aver- 
age standard deviation of the flux at two windows: 1450 A - 
1470 A, and 1710 A - 1730 A (rest frame) 1 . The rest frame 
wavelengths were calculated using redshifts determined from 
the peak of the [O III] A5007 line (kindly provided by T. 
Boroson, private communication). The wavelength depen- 
dence of the flux density uncertainty (of the HST spectra) 
was assumed to scale with the square root of the ratio of 
the flux density to continuum flux density. The IUE spectra 
were obtained with the flux density uncertainty from the 
IUE archives. 

The spectra were corrected for Galactic reddening us- 
ing the E(B— V) values from Schlegel, Finkbeiner & Davis 
(1998, as listed in the NA SA /IPAC Extr agalactic Database), 
and the reddening law of ISeato 3 (Il979l) . A local power-law 
continuum was fit to each spectrum between ~ 1470 A and 



1 Lack of data or strong systematic features necessitated the fol- 
lowing changes to one or both windows; PG 0003+158: 1465 A - 
1485 A, PG 1114+445: 1465 A - 1485 A, and 1730 A - 1750 A; 
PG 1211+143, PG 1404+226 and PG 1440+356: only the 1710 A- 
1730 A window used; PG 1415+451: used only 1700 A - 1720 A; 
PG 1543+489: 1600 A - 1620 A; PG 1612+261: 1470 A - 1490 A; 
PG 2304+042: 1680 A - 1700 A. 



~ 1620 A, and the C IV line emission was fit as a sum of 
three Gaussians, using the procedure described in Laor et 
al. (1994, §3 and the appendix there). Wavelength regions 
suspected to be affected by intrinsic or Galactic absorption 
were excluded from the fit (Laor & Brandt 2002). The pur- 
pose of the line fit is not to decompose the line to possible 
components, but rather to obtain a smooth realization of 
the line profile, which is likely to yield more accurate values 
for the line width and EW. 

Table 1 presents our measured rest-frame EW of the 
best-fitting models to the C IV profile, together with the es- 
timated errors. The errors were estimated by repeating the 
model fits with the power-law continuum displaced upward 
and downward by la. As expected, the typical errors asso- 
ciated with IUE spectra are significantly higher than those 
of the HST spectra. One should note that larger systematic 
errors could be associated with the placement of the con- 
tinuum windows. In particular, the presence of very broad 
weak wings is difficult to detect, although they may have 
a non-negligible EW. Table 1 also lists vL v calculated at 
3000 A for each object, and the estimated L/L^d- 



3 THE CORRELATION ANALYSIS 

Figure 1 (top panel) presents the Baldwin effect for the 81 
BQS quasars measured above. The C IV EW are taken from 
Table 1, and for the luminosity we use vL v at 3000 A (see 
§2). To verify that our sample is not biased in any way, we 
also present in Fig. 1 data from two other stud ies, 454 Large 
Bright Quasar Survey objects 2 analyzed by iForster et alJ 
fcOOll). and 125 pre-COSTAR AGN observations analyzed 
bv iKuraszkiewicz et al] i2002f) (which overlap some of our 
objects). These data are available in public electronic form. 
The luminosities of these objects were derived from the spec- 
ified flux densities using the same cosmology assumed for the 
BQS quasars. The trend and scatter of the Baldwin effect 
for the BQS quasars appears to be ver y similar to those dis- 
played by the lForster et all i200ll) and IKuraszkiewicz et al] 
(2002) samples. The Spearman rank-order correlation coef- 
ficient for the BQS Baldwin effect is r s = —0.154. The flat- 
tening of the Baldwin relation at low luminosity is clearly 
seen in Fig. 1 (see discussion in §4.2 of Osmer & Shields 
1999). The non-linearity of the Baldwin effect (on log-log 
scale) suggests that the Spearman rank-order correlation 
coefficient, which tests for a monotonic relation, is a more 
suitable statistical test here, compared to the Pearson cor- 
relation coefficient, which tests the significance of a linear 
relation. 

Table 2 presents the results of a correlation analysis 
of the C IV EW reported in Table 1 with all the optical 
emission line parameters from BG92 (table 2 there), sup- 
plemented by ot ox from Brandt, Laor & Wills (2000), and 
L/LEdd and Mbh as estimated from the optical luminos- 
ity and the H/3 FWHM (equation 3 in Laor 1998). For the 
sake of brevity we present only the most significant correla- 
tions (those where Pr < 5 x 10~ 4 ). To explore the effect of 
the generally lower S/N IUE spectra on the strength of the 

2 We have i gnored two upper limits in the tabulation of 
IForster et alJ 1200 lfl . and several obvious typographical flux er- 
rors, in their total list of 488 objects. 
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Table 1. The Civ EW, vL u , and L/L Edd of the 81 BQS quasars. 
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"■ In units of log(^Z/„/10 44 ), where uL v is measured in erg s — 1 at rest frame 3000A. 
c Log L/L Edd . 



correlations, we also present in Table 2 the correlations ob- 
tained based only on the 46 objects with the generally higher 
S/N HST spectra. The correlations for the HST sample are 
similar, which suggests that the C IV EW measurement un- 
certainty is not a major source error. 

The strongest correlations of the C IV EW are with a ox , 
and with parameters related to the strength of the [O ill] 
and Fell lines, and with the H/3 FWHM, which are part of 
the BG92 EV1 correlations. These correlations confirm ear- 
lier results on the inverse relation of the C IV EW and the 
Fen strength (Marziani et al. 1996; Wang et al. 1996), the 
Civ EW EV1 correlations found by Wills et al. (1999, table 
1 there) for a smaller sample of 22 BQS quasars (part of 
our sample), and the earlier results on the C IV EW correla- 
tion with ot ox (see §1). The correlation with a ox is consistent 
with the suggestions that the Baldwin effect is driven by a 
softening of the ionizing continuum (which decreases a ox ), 
with increasing luminosity (Korista et al. 1998) . The inverse 
correlation of the Civ EW and the Fen strength may be 
due to optical depth effects. A large optical depth in the Fe n 
emitting region strengthens the optical Fe n emission by con- 
verting UV Fe 11 emission to optical Fe 11 emission (Netzer & 
Wills 1983, Shang et al. 2003). In contrast, a large optical 
depth in the ClV emitting region results in collisional sup- 
presses of this line (Ferland et al. 1992) . Thus, the observed 
inverse correlation could be produced if the optical depths 
of the Fe n and C IV emitting regions are related. 



Objects with a ox < -2, or 'Soft X-Ray Weak' AGN, 
generally show intrinsic broad C IV absorption (Brandt et 
al. 2000), whose strength is correlated with luminosity (Laor 
& Brandt 2002). Incomplete correction for this absorption 
may induce the observed trend of decreasing C IV EW with 
decreasing a ox . To avoid absorption biases, we repeated the 
C IV EW vs. a ox correlation excluding the 12 objects with 
a ox < —1.75 where there is potentially significant absorp- 
tion. This gave essentially the same result (rs = 0.498, vs. 
rs = 0.525 for the complete sample), indicating there is no 
significant absorption bias. 

The interesting new result in Table 2 (Fig. 1, middle 
panel) is the strong correlation of the C IV EW with the 
L/L Edd indicator (i.e. L 1/2 (H/3 FWHM) -2 ), where r s = 
-0.581 (null probability of Pr = 1.3 x 10" 8 ). This is not 
a simple consequence of the separate correlations of the 
Civ EW with L and with the H/3 FWHM, as for exam- 
ple, the correlation with Af BH (i.e. L 1/2 (R/3 FWHM) 2 ) is in- 
significant (rs = 0.21, Pr = 0.05). However, does the combi- 
nation of L and Hb FWHM in the form L 1/2 (H/3 FWHM)" 2 
produces the tightest correlation? To answer that we 
searched for the largest rs in a correlation of log C IV EW 
with the linear combination log L+a log(H/3 FWHM), where 
a is a free parameter. We found a maximum correlation of 
rs = —0.588 at a = —5.3, close to the combination of L 
and H/3 FWHM which provides L/L Edd (a = -4). Thus, 
L/L Edd may be the primary parameter which drives the 
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Baldwin effect for C IV. This implies that the large scatter 
in the Baldwin effect is produced by the range of L/LEdd at 
a given L, as suggested by Shang et al. (2003). However, our 
results disagree with the suggestion of Shang et al. that the 
ClV EW is primarily dependent on L, rather than L/LEdd- 

Is L/Z/Edd the only parameter which controls the 
C IV E W? Since the EV1 correlations also appear to be 
driven by L/LEdd (BG92; Boroson 2002), their correlation 
with the C IV EW could be interpreted as secondary corre- 
lations. This can be tested with a partial correlation analy- 
sis, which yields an average rs = 0.42 for the correlation of 
C IV EW with a ox and the EV1 parameters, keeping L/Z/Edd 
fixed, and also for the correlation of Civ EW with L/Z/Edd, 
keeping EV1 and a ox fixed. Although these correlations are 
weaker than the original ones (Table 2), they are still signifi- 
cant (Pr ~ 10 -4 ). The significance of the partial correlations 
suggests there is a true scatter in the BLR properties at a 
fixed L / LEdd , though it may also be induced by the inaccu- 
racy of our L/Lsdd indicator. 

The significance of the partial correlations implies that 
the scatter in the Civ EW vs. L/LEdd relation can be re- 
duced by including a third parameter (in addition to L 
and the H/3 FWHM, which form L/L E dd). The param- 
eter which most significantly improves this correlation is 
the [Om] A5007 EW (Fig. I, lower panel), which yields 
r s = -0.722 (using a ox instead of [O ill] A5007 EW as a 
third parameter yields r s — —0.667, and using the Fell/H/3 
flux ratio yields r 3 = —0.674). Similarly, the most signifi- 
cant fourth parameter is cv ox (i.e. correlating log C IV EW 
vs. a linear combination of log L/LEdd, log [O ill] A5007 EW, 
and a ox ), which yields a small improvement to r s — —0.753. 
Finally, the inclusion of the Fell/H/3 flux ratio as a fifth 
parameter yields a slight improvement to r s — —0.767. 

Some of the scatter in the C IV EW correlations may 
be induced by a spread in the BLR covering factors among 
AGN. A plausible way to overcome such a spread is to use 
the Civ EW/FL3 EW ratio in the above correlations, in- 
stead of the C IV EW. This yields significantly lower corre- 
lations (e.g. the correlations with L/LEdd and a ox go down 
to -0.391 and 0.391 from -0.581 and 0.525, respectively). 
This may suggest that the H/3 EW is not a good indicator of 
the BLR covering factor, but is rather modulated by other 
BLR properties which do not modulate the C IV EW. 



4 CONCLUSIONS 

The purpose of this paper is to obtain some indications for 
the physical parameters which drive the Baldwin effect in 
C IV through correlation analysis. We use archival HST and 
IUE C IV spectra of sufficient quality available for 81 of the 
87 BQS quasars, together with the optical emission line pa- 
rameters from BG92, and a ox from Brandt et al. (2000). 

We find that a major source of scatter in the Baldwin ef- 
fect is the H/9 FWHM. Its inclusion as a second parameter to 
L, in the form L 1/2 x(H/3 FWHM) -2 , which is proportional 
to L/LEdd, leads to a significant increase in the correlation 
strength (from r 3 — —0.154 to r s = —0.581). This indicates 
that the Baldwin effect is a secondary relation, which is in- 
duced by the stronger relation of the C IV EW and L/L E dd, 
and the tendency of high L AGN to have a higher L/Lsdd- 
This also explains why NLSIs, which are high L/LEdd AGN 
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Figure 1. The main Civ EW correlations discussed in this pa- 
per. The Spearman rank-order correlation coefficient (rs) and 
the null probability (Pr) are indicated at each panel. Top panel: 
The Baldwin effect for the 81 BQS quasars , toget her with the 
iForster et alj 1200 ll) and lKuraszMewicz et alj 120021) samples. All 
samples are similarly distributed. Middle panel: The correlation 
of the Civ EW with L/L E dd- Note the significant increase in rg. 
Lower panel: The correlation with the addition of [O III] EW as a 
third parameter, which further increases rg. 
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Table 2. The main Civ EW correlations. 



Variable Name" 


h 

r a ° 






i/L„(3000A) 


-0.154 


1.71x10" 


-01 




—0.018 


9.08x10" 


-01 




-0.581 


1.31x10" 


-08 




—0.642 


1.53x10" 


-06 


Ql ox . . 


0.525 


4.87x10" 


-07 




0.463 


1.18x10" 


-03 


[Om] A5007 EW 


0.624 


4.71x10" 


-10 




0.708 


3.67x10" 


-08 


Fen EW 


-0.518 


7.49x10" 


-07 




-0.536 


1.24x10" 


-04 


H/3 FWHM 


0.427 


7.03x10" 


-0B 




0.510 


2.92x10" 


-04 


R [O III] A5007 peak height 


0.624 


4.78x10" 


-1(1 




0.647 


1.20x10" 


-06 


R Fe II EW 


-0.626 


4.02x10" 


-10 




-0.698 


6.94x10" 


-08 


R [O in] A5007 EW 


0.471 


9.23x10" 


-06 




0.494 


4.89x10" 


-04 



a The prefix 'R' indicates the ratio of the variable to that of H/3 
(BG92), e.g. R Fell EW corresponds to Fell EW/H/3 EW. 
" The top value for each variable is obtained for the complete 
sample of 81 objects, the bottom value corresponds to the sub- 
sample of 46 objects with HST spectra. 



at a low L, have much lower C IV EW than expected for 
their L (as found by Wilkes et al. 1999). 

We do not know the physical mechanism responsible for 
the reduction in the C IV EW with increasing L/L-Edd- How- 
ever, we find that the scatter in the Baldwin effect can be 
further reduced by including either the [O ill] EW, the rela- 
tive Fell strength, or a ox , as a third parameter. The Fell 
strength may be an indicator of the BLR optical depth, 
and a x may be an indicator of the ionizing continuum 
shape, which are known theoretically to affect the C IV EW. 
The physical mechanism implied by the correlation with the 
[O ill] EW remains to be understood. At least some of the 
remaining scatter can be produced by a non-isotropic, or a 
time variable continuum source. 

The results presented here can be tested with larger 
samples of AGN. In particular, it will be interesting to ex- 
plore the extension to high z AGN, where the H/3 region 
is observable in the IR (e.g. Mcintosh et al. 1999; Yuan & 
Wills, 2003; Shemmer et al. 2004). In addition, different lines 
show different slopes for the Baldwin effect, and it will be 
interesting to explore what are the primary physical param- 
eters which drive other emission lines. 
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